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Abstract. Inspiral signals from binary black holes, in particular those with 
masses in the range IOMq $ M ^ IOOOMq , may last for only a few cycles within 
a detector's most sensitive frequency band. The spectrum of a square-windowed 
time-domain signal could contain unwanted power that can cause problems in 
gravitational wave data analysis, particularly when the waveforms are of short 
duration. There may be leakage of power into frequency bins where no such power 
is expected, causing an excess of false alarms. Wo present a method of tapering 
the time-domain waveforms that significantly reduces unwanted leakage of power, 
leading to a spectrum that agrees very well with that of a long duration signal. 
Our tapered window also decreases the false alarms caused by instrumental and 
environmental transients that arc picked up by templates with spurious signal 
power. The suppression of background is an important goal in noise-dominated 
searches and can lead to an improvement in the detection efficiency of the search 
algorithms. 
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1. Introduction 

Interfcrometric gravitational wave detectors are now operating at sensitivity levels 
at which one can expect to detect inspirals from compact binary coalescences at the 
rate of one per 5 years (optimistic rate) to one per 5000 years (pessimistic rate). 
When upgraded to advanced detectors, these might be as large as 400 per year to 
one per 2.5 yrs [1]. Even so, most of the inspiral signals are not likely to stand 
above the noise background. A variety of techiques to enhance signal visibility and 
reject false alarms are currently being used in gravitational wave searches. Examples 
include matched filtering for signals of known phase evolution [2] , wavelet transforms 
for transient signals of unknown shape [31 |3], coherent search methods for burst 
signals [5], etc. Moreover, vetoes based on the expected signal evolution [B] and 
instrumental and environmental monitors [7] have been developed over the past decade 
to improve detection probability and mitigate false alarms. Detecting a signal buried 
in non-stationary noise is a challenging problem as some types of non-stationary noise 
artefacts can partially mimic the signal. 
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Many of these techniques involve the computation of a correlation integral in 
which band-passed data are multiplied by the frequency-domain (FD) model waveform 
or the discrete Fourier transform (DFT) of the time-domain (TD) signal (see, for 
example, Here we consider a matched filtering search for inspiral signals where 
the DFT of a TD waveform is used to construct the correlation. A problem that has 
not been adequately addressed (see, however, Ref. [9]) in this context is the effect of 
the window that is used in chopping a TD signal before computing its DFT. 

Inevitably, all signal analysis algorithms use, implicitly or explicitly, some form 
of window function. An inspiral waveform sampled from a time when the signal's 
instantaneous frequency enters a detector's sensitive band until the time when it 
reaches the frequency of last stable orbit (FLSO) implicitly makes use of a square 
window. Signal analysis literature is full of examples of artefacts caused by the use of 
such window functions. For instance: leakage of power from the main frequency bin 
where the signal is expected to lie into neighbouring bins, loss of frequency resolution 
and corruption of parameter estimation [TU] . In this paper we explore the problems 
caused by using a square window and suggest an alternative that cures some of the 
problems. 

There is no unique, or favoured, windowing method. One is often guided by 
the requirements of a particular analysis at hand. In our case, a square window is 
especially bad since the leakage of power outside the frequency range of interest can 
lead to increased false alarm rate and poorer estimation of parameters. One reason for 
increased false alarm rate could be that the noise glitches in the detector look more 
like the untapered waveform and less like a tapered one. We have explored the effect 
of a smoother window function, presented in Section [21 which has a far steeper fall-off 
of power outside the frequency range of interest. Use of this window has cured several 
problems we had with a square window. In Section |4] we will discuss how tapering 
helps in a more reliable signal spectral estimation and hence a proper determination 
of the expected signal-to-noise ratio. Spectral contamination is worse for larger mass 
black hole binaries as they are in the detector's sensitive band for a shorter time and 
the window function can only extend over a short time. It is for such signals that 
our tapered window offers the most improvement. In Section [S] we will discuss how 
the rate of triggers from a matched filtered search can vary depending on the kind of 
window function used. We shall briefly mention in Section |6] what effect our window 
function has on parameter estimation, giving the conclusions of our study in Section[71 

2. Window functions and their temporal and spectral characteristics 

Let h{t) denote a continuous differentiable function, for example a gravitational wave 
signal emitted by a coalescing compact binary, and let H{f) denote the Fourier 
transform (FT) of h{t) defined by 



In reality the signal really does not last for an infinite time. The FT of a signal of 
finite duration lasting, say, from —T/2 to r/2, can be represented either by using the 
limits of the integral to go from —T/2 to T/2 or by using a window function. The 
latter is preferred so as to preserve the definition of the FT. 

A window function is a function that has either a finite support or falls off 
sufficiently rapidly as t — > ±oo. Two simple windows that have finite support are 




(1) 
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the square window st(0 defined by 



1, -^<t<^ 



- 1 O! otherwise: 
and the triangular window 6t(0 defined by 

Neither the square nor the triangular window are differentiable everywhere. As a 
result, they are not functions of finite bandwidth. In other words, their FTs, 5t(/) 
and -Bt(/), do not have finite support in the FD: \S{f)\ > for — oo < / < oo. In 
the case of a square window the FT S{f) is a sine- function, |5t(/)| = rsinc(7r/r), 
which is peaked at / = 0, with a width tt/T and falls off as as / — ?> ±oo. The lack 
of finite support in the Fourier domain could sometimes cause problems, especially 
when the width of the window in TD is too small. For functions that have infinite 
bandwidth the sampling theorem does not hold but this is not a serious drawback if 
the FT falls off sufficiently fast above the Nyquist frequency. However, there could be 
other issues when the window leads to leakage of power outside a region of interest as 
we shall see below. 



2.1. The Planck-taper window function 

A signal h{t) with the window WTit) applied to it, in other words the windowed signal 
hw{t), is defined by 

h^{t) = h{t)wT{t). (4) 

The convolution theorem states that the FT of the product of two functions h{t) and 
WTit) is the convolution of individual FTs: 



Hu,if) 



h{t)wT{t) exp(27ri/t)dt 



= H{f)*WT{f) 



HinWrif - f'W. 



(5) 



(6) 



We can now see why a window whose power in the FD does not fall off sufficiently 
rapidily might be problematic. The convolution integral will have contributions from 
all frequencies. Suppose we are interested in matched filtering the data with an inspiral 
signal from a compact coalescing binary whose instantaneous frequency varies from fa 
at time ta to at time t},. One would normally achieve this by using a square window 
ST{t) that is centered at [ta + tb)/2 with width T = tt — ta- However, we can see from 
(j6|) that the convolution integral will have contributions from outside the frequency 
range of interest. 

To circumvent this problem we propose to use a window function that falls off 
rapidly outside the frequency range of interest. Inspired by the tapering function used 
in Damour et al we define a new window function a(t) by 
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Figure 1. The Planck-tapcr window in the TD (left), for three different ehoices of 
the parameter e = 0.01, 0.033, 0.1, and their power spectra (right). For reference 
we have included the square window with the same effective width as the Planck- 
taper window. 



Here T is the width of the window and e is the fraction of the window width over 
which the window function smoothly rises from at i = ti to 1 at t = i2 or falls from 
1 at t = to at i = t^. We shall call ait) the Planck-taper window as the basic 
functional form is that of the Planck distribution. The motivation for choosing this 
window function is to reduce the leakage of power in the FD but at the same time 
not to lose too much of the length of the signal in the TD. The choice of e will affect 
both aspects significantly. In figure [T] we have shown the window function for several 
choices of the parameter e = 0.01,0.033,0.1.. At lower frequencies the spectrum of 
the Planck-taper window falls off at the same rate (i.e., 1//) as a square window. But 
beyond a certain frequency fo {eT)~^, the spectrum falls off far faster. 

A key feature in our use of the Planck-taper window is the automated and 
waveform-dependent adjustment of e as discussed in section [2 . 2 1 below. 



2.2. Implementation of the window 

We discretise (O by replacing t, ti, ^2, ^4 with the array indices j, ji, j2, js, j4- In 
this notation the parameter epsilon is approximated by e ~ (j2 — ji)/N, where N 
is the number of data points in the waveform. The start and end of the waveform 
are denoted by ji and j4, respectively. The values of j2 and have to be chosen 
judiciously to avoid leakage of power. We choose j2 and ja to be the array index 
corresponding to the second stationary point after ji and before ji (see figure [5]). 
Applying the transition stage of a from a crest /trough ensures that the window does 
not have a sudden impact on the behaviour of the waveform. The first stationary 
point would not be an appropriate choice as it may occur within only a few array 
points of ji or j4, causing e to be too small. One could choose the 3rd, 4th or 5th, 
but using such later maxima would reduce the genuine power of the waveform more 
than what might be acceptable. 
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Figure 2. The window function has 
been applied to the start of a cosine 
wave (top curve) using two methods. 
In the first case it is applied from 
J = 1 up to an arbitrary choice of 
j = 20 (middle) , whereas in the second 
case it is applied up to the second 
maximum at j = 100 (bottom). The 
lighter coloured parts of the middle 
and bottom curves (to the left of the 
black vertical lines) show where the 
taper has been applied. 



2.3. Comparison with other windows 

We do not compare the performance of Planck-taper with other commonly used 
windows, e.g., Bartlet, Hann or Welch. Such windows transition between and 1 over 
j = 1, . . . ,N/2, where the window is of length N, producing significant differences 
between h{t) and hwit) in (|4]). The power is therefore supressed at the beginning and 
end of h{t). This is acceptable when computing the power spectral density (PSD) of 
a data segment, but would cause a problem if applied to a template waveform as the 
phase (frequency) and amplitude of h{t) are both instantaneous functions of t, with 
the most power at the end of the waveform. 

Windows with properties similar to Planck-taper, such as having a central flat 
region, do exist. For example, the Tukey window |12j . which has been used in 
gravitational- wave data analysis recently [13], may offer a good comparison. However, 
a key feature in our study of the Planck-taper window is the waveform-dependent 
adjustment of j2 and .73. Whilst this automation could be considered separately from 
the Planck-taper window and used on other windows defined by the points ji,2,3,4, we 
have not done so here. Given the shared features of the Tukey window with Planck- 
Taper one might expect similar results. 

3. Effect of the window^ function on the signal spectrum 

In this section we will examine the power spectrum of the waveform of a coalescing 
binary emitted during the inspiral phase. The waveforms arc modelled using the post- 
Newtonian (PN) approximation. However, even within the PN approximation, there 
are several different ways in which one might construct the waveform |14| I15j . Two 
such models widely used in the search for compact binary coalescences are Taylor 
T3 (TT3) and the stationary-phase approximation (SPA). TT3 is a TD signal model 
in which the amplitude and phase of the signal are both explicit functions of time. In 
the so-called restricted PN approximation the signal consists of the dominant harmonic 
at twice the orbital frequency, but not higher order PN corrections consisting of other 
harmonics, and the phase is a PN expansion that is currently known to 0{v'^) in the 
expansion parameter v - the relative velocity of the two stars. The SPA is the Fourier 
transform of the TT3 model obtained by using the stationary phase approximation 
to the Fourier integral [16]. A template belonging to the TT3 model is defined for 
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times when the gravitational wave frequency is within the detector's sensitivity band 
until it reaches FLSO. This means one is in effect multiplying a square window with 
a continuous function. 

Figure |3] shows the signal-to-noise ratio (SNR) integrand of the SPA, computed 
using the initial Laser Interferometer Gravitational-wave Observatory (LIGO) design 
PSD [14]. The inspiral waveform is defined from a lower cutoff frequency of 35 Hz 
up to its FLSO, for 20 Mg and 80 Mg equal-mass binaries. The DFT of the TT3, 
generated between the same frequencies, with a square window (or rather no window), 
labelled H5 , and with the Planck-taper window, labelled H,^ , are also plotted. Where 
the Planck-taper window is used the excess power, that above FLSO, decreases rapidly 
and the spectrum is closer to that of the SPA. 



20 Mb 80 M( 




10^ 10^ 10^ 10^ 

Frequency /Hz Frequency /Hz 



Figure 3. The plots show the SNR integrand, where the waveform is generated 
from a frequency of 35 Hz to the FLSO of the source, computed using the initial 
LIGO sensitivity for sources of total mass 20 Mq and 80Mq . In the case of the 
Planck-taper window the SNR integrand falls off far faster than in the case of the 
square window or the TT3 model. 



4. Effect of the window function on the estimation of the signal-to-noise 
ratio 

Gravitational wave searches for known signals, such as those emitted from a compact 
binary coalescence (CBC) [171 [18], rely upon signal models for two primary reasons. 
Firstly, they are used as templates to matched filter the data. Secondly, they are 
injected into the data as simulated signals to estimate the efficiency of the detector to 
detect such signals. If the signal/template models are generated in the TD then they 
must undergo a DFT if the data are analysed in the FD as is the case for the current 
LIGO matched filter code. 

The expectation value for the SNR of a signal in stationary Gaussian noise, when 
the signal and template match exactly, is given by 

where Sh{f) is the detector noise PSD and f lower is the detector lower cutoff frequency 
chosen so that the contribution to the SNR integral from frequencies / < fiower is 
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negligible. The second of the expressions on the RHS is a discretised evaluation of the 
SNR which is often used in numerical calculations. Here H^, k ~ 0, . . . , N/2, is the 
DFT of the signal defined for positive frequencies and S^k is the discretised PSD. 

The amplitude of an inspiral signal increases with the total mass of the system; 
conversely, the FLSO of the signal is inversely proportional to the total mass. 
Therefore, as the total mass of a system increases, the amplitude of the signal and the 
FLSO will have opposing effects. For lower mass systems, the increasing amplitude 
causes the SNR to increase as a function of the total mass. However, for higher mass 
systems, the reduction in the FLSO causes the signal to have less power in band. As 
a result, the SNR will decrease as a function of the total mass. The relatively low 
FLSO of the higher mass templates, coupled with their short duration, lead them to 
be particularly susceptible to artefacts of spectral leakage in the DFT. 

Figure 2] shows the SNR for TT3 inspiral waveforms that are 2PN in amplitude 
and phase, plotted as a function of the total mass for two choices of the window 
function: the dashed curve corresponds to the square window and the solid curve to 
the Planck-taper window. All other parameters are the same in both cases. When 
the Planck-taper window is used, the curve exhibits the expected behaviour, whereas 
in the case of a square window , the SNR curve is 'jagged' which is unexpected given 
that stationary Gaussian noise was used in the estimation of the SNR. This behaviour 
is most likely explained by the excess power from the DFT of the waveform. 
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Figure 4. The SNR vs the total mass 
of the source for signals correspond- 
ing to compact binary systems directly 
overhead a detector of initial LIGO de- 
sign PSD. We plot the SNR obtained 
using the DFT of TD waveforms with 
a square window (dashed curve) and 
with the Planck-taper window (solid 
curve) . Here the systems are overhead 
the detector at an effective distance of 
65Mpc, using a fixed mass ratio of 5 : 1 
and a fixed inclination angle of 45°. 



It should be noted that integrating to FLSO rather than Nyquist in Eq. ([5]), is not 
considered appropriate here. Firstly, the higher harmonics in the amplitude corrected 
waveforms contain power above FLSO, (which becomes more signifcant for high mass 
systems). Secondly, cutting off the integration at FLSO is essentially the application 
of a square window to the template waveform in the frequency domain. This will 
lead to leakage of power in the time domain which is not a desirable feature. The 
problem of using a square-windowed TD template as our matched filter is not that 
there is power above FLSO; it is that the excess power in this region, present due to 
windowing, but not present in a genuine signal will lead to unnecessary false alarms 
in a search. 
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5. Effect of window functions on Trigger Rates 

To assess the effect that tapering of templates has on trigger rates, we have apphed the 
LIGO Scientific Collaboration (LSC) CBC pipeline [HI EQl [U UHl EJ to data taken 
during the fourth science run (S4) of the LIGO, which took place from February 22 
- March 23, 2005. The basic topology of the pipeline is similar to that used in many 
previous searches [20l [171 [18] , and consists of the following main steps: 

• The template bank is chosen such that the loss of SNR due to having a finite 
number of templates is no more than 3% for any signal belonging to a given 
family of waveforms [HI [2] • 

• Matched filter the data with the generated templates. A trigger is generated at 
times when the SNR is larger than a given threshold. The output of this stage is 
a list of first-stage single- detector triggers. 

• Check for coincident events between different detectors. For an event to be 
deemed coincident, the parameters seen in at least two detectors (for instance, 
the masses of the system, the time of coalescence, . . . ) should agree to within a 
certain tolerance [23]. The output of this stage is a list of first-stage coincident 
triggers. 

• Re-filter the data using only templates associated with coincident triggers. This 
time, the triggers are subjected to further signal-based vetoes, some of which are 
computationally costly, such as the chi-squared veto [6]. This produces a list of 
second-stage single- detector triggers. 

• Check for coincident events between detectors using the second-stage single- 
detector triggers. This produces a list of second-stage coincident triggers. 

In this study the data were filtered using the effective one-body (EOB) templates 
[24l [25l IT4] , tuned to recent results in numerical relativity |26l [27] , with a total mass 
in the range 25 — lOGM©. This choice agrees with the templates used to search for 
signals from high-mass CBCs in data from LIGO's fifth science run (S5). Because the 
EOB waveforms used as templates contain the inspiral, merger and ringdown phases, 
there was no need to taper the end of the waveform. Therefore, in this case, the taper 
specified in ([7]) was only applied to the start of the waveform. Although this may 
reduce the effect the taper has in comparison to tapering both ends of an inspiral- 
only template, it is of more interest to evaluate the performance in a realistic search 
case. It should be noted that the tapering window is explicitly applied to the template 
waveform where the length of the waveform is less than the length of the data segment 
that is matched filtered. We do not apply any window to the data segment. 

Figure[5]shows the number of triggers as a function of total mass with and without 
tapering. It can be seen that the number of triggers is generally higher when the 
templates are not tapered. The only exception seems to be the lowest mass bin in the 
second-stage coincident triggers, where the opposite is true. However, the difference in 
the number of triggers in this bin is not large, and is likely just a statistical anomaly. 
For first-stage single-detector triggers, the number of triggers using tapered templates 
is 84% of that obtained using un-tapercd templates. The number of second-stage 
coincident triggers when using tapered templates is 71% of that obtained for un- 
tapered templates. The difference in trigger rates is more significant at higher masses. 
This is because the template waveforms for these systems terminate at a frequency 
within or below the most sensitive frequency band of the detector, making any leakage 
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of power to higher frequencies more significant (cf Figure [31 left most panel). The 
reduced trigger rate indicates that applying the taper function to the templates could 
aid in reducing the false alarm rate in a search for high mass CBCs. 
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Figure 5. Left: Number of triggers reeovered by matclied filtering the S4 data 
witli and without tapering applied to the templates. Right : Number of triggers at 
the second stage of the analysis pipeline, after consistency checks and coincidence 
tests 1231 in the time-of-arrival and masses of the component stars have been 
applied. 



6. EfTect of windowing on detection efficiency and parameter estimation 

The same data used in section [5] were re-analysed, but with simulated gravitational 
wave signals (injections) added. The injections were of the same family as the 
templates used in section [5j This allowed us to compare the detection efficiencies 
and accuracy of parameter estimation using tapered vs. untapered templates. We 
looked at the error in recovered chirp mass and arrival time at both single detector 
first stage triggers and coincident second stage triggers, but found negligible difference 
between the two cases. 

We did not explicitly measure the detection efficiency as a function of distance, 
but found the number of injections recovered to be nearly identical in the two cases, 
with less than 1% fewer injections found when using tapered templates. Given the vast 
reduction in the trigger rates shown in Section [SJ this indicates that an improvement 
in detection efficiency can be expected when using tapered templates. 

The above studies were performed first with tapering applied to the injections 
and then repeated without - the difference between the results was negligible. 

7. Conclusion 

We have developed a tapering method that leads to a spectrum for TD waveforms 
that more closely matches their FD analogs, containing significantly less power at 
unexpected frequencies when compared with the use of a square window. This is 
achcivcd by automating the implementation of the window. 

I We have seen some evidence of improvements in parameter estimation for the ambiguity function 
of high mass inspiral-only waveforms, but this is outside the context of a gravitational wave search. 
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If tapering is applied to templates in a gravitational wave search the trigger rates 
are reduced, especially for high mass templates, without any significant change in 
detection efficiency. In a search, foreground triggers can be ranked by their probability 
of occurring as a background trigger; thus if background triggers are reduced, a given 
foreground trigger may appear more significant. Another benefit of reduced trigger 
rates is that the computational cost of a search will decrease. We have demonstrated 
that the windowing method would be beneficial when used in a high mass search. 

The tapering method could also be useful in low latency data analysis techniques 
where TD templates are divided into sub-templates of different frequency ranges, and 
matched filtered individually [28]. The relative shortness of some templates in the 
higher frequency bands potentially compounds the problem of using a square window, 
and tapering the templates may go some way to alleviating this issue. 
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